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SUMMARY 

A coaiputatlonal Inverse procedure for transonic airfoils In which shapes 
are deteroilned supporting prescribed pressure dlsuributlons Is presented. 

The aiethod uses the soiall disturbance equation and a consistent analysis- 
design differencing procedure at the airfoil surface. This avoids the 
Intenaedlate analysis— design— analysis Iterations. The effect of any openness 
at the trailing edge Is taken Into account by adding an effective source term 
In the far field. The final results from a systematic expansion procedure 
which models the far field for solid. Ideal slotted and free Jet tunnel walls 
are presented along with some design results for the associated boundary 
conditions and those for a free flight. 

INTRODUCTION 

Computational design or Inverse procedures for transonic airfoils In 
which shapes are determined supporting prescribed pressure distributions have 
been In use since the early work of Nleuwlandl** which employed hodograph 
methods to calculate shock— free supercritical flow about a family of quasi— 
elliptical airfoils. Later Garabedlan and Korn^ developed a more general 
hodograph procedure to design highly cambered shock— free airfoils. In spite 
of their usefulness, hodograph procedures for design purposes have several 
disadvantages. They require too many Input parameters, are restricted to 
shock-free solutions, and are not easily extendable to design of three-dimen- 
sional wings. Steger and Kllneberg^ treated the problem within a small- 
dlsturbance framework solving the continuity and vortlclty equation at 
Interior points. To Insure consistency between the analysis and the design 
formulation, they applied appropriate discretization procedures to the 
vortlclty equation at the airfoil grid points. However, the first-order 
system with velocity components as dependent variables produces a difficulty 
In the treatment of singularities at the airfoil nose and trailing edge. The 
effect of nose and trailing edge singularities could be greatly reduced by 
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f foj^latlon involving the velocity potential. Tranen^ employed 
the full potential equation to remedy the deficiency Inherent in the small 
disturbance formulation at the leading and trailing edges. To overcome the 
inaccuracies and inconsistencies in his formulation associated with the 
discretization procedures at the boundary, iterations must be employed between 
direct and Inverse solvers. Also, in the full potential formulation, boundary 
conditions are to be applied at the exact airfoil surface. Since the airfoil 
surface is unknown in the design problem, errors propagate due to application 
of boundary conditions at some assumed airfoil surface. Carlson^ used ghost 

accurate methods to handle airfoil boundary points. 
®“Ploylng the circle plane as in Tranen’s procedure, Carlson used 

Je^M thif w T” •'i® procedure is not consistent in the 

^ discretizations used for Cp in the analysis and those in the 

a ®®“® former uses a central differenced 

♦x Cp calculations while the latter employs a special backward differenced 
ghost points. As a result, perfect agreement between his analysis 
calculations is not to be expected, especially near the shock. 

Mso, the problems associated with open trailing edges are not addressed in 
his work. The trailing edge is made to close by altering the nose shape. 

*‘^® intermediate calculations, large open trailing edges occur. 

®'® ‘!"®®tio“ahla since the effect of openness is not 
Included in the far field for the intermediate solutions. 

in thl® pap®r, a small-disturbance model employing the velocity potential 
as the dependent variable is used for the implementation of the design 
algorithm. This procedure simplifies the treatment of boundarj' conditions 

f "fj “«PPing8 that arise in a full potential equation 
formulation. A mixed boundary value problem is solved in which Neumann data 

^®” P®^®®nt oi the chord length where the assumed 
flirfM?* ®”‘^ Dirichlet conditions are prescribed on the rest of the 

airfoil where the pressure is to be modified. One important thrust of the 
airfnii developing a consistent discretization procedure for the 

points. If the converged Cp output from the analysis is not 
altered, then the design mode recovers the same airfoil shape without any 

In^Carlion’a^ In the airfoil slope at the shock wave, overcoming a deficiency 
in Carlson s work. Another significant feature of the numerical implementa- 
tion not considered by the previous investigators is the effect of an open 

with Jhf field. In the present work, this is accounted for 

addition of the necessary source terms in the far field. Desisn of 
thick trailing edge airfoils is of interest in invlscid flow to achieve a 
r®asonable trailing edge thickness after accounting for the viscous displace- 
ment thickness. Some amount of trailing edge thickness is required from a 
structural stability point of view. Figure 1 schematically explains the 
design philosophy followed in this paper. The top of Figure 1 shows a conven- 
tional airfoil at transonic speed producing a shock on the upper surface, 
pec fying a shockless pressure distribution on the upper surface would 
flatten the upper surface of the conventional airfoil, thereby producing an 
openness at the trailing edge. This is shown in the middle of Figure 1. The 
amount of trailing edge openness can be reduced by specifying a lower surface 
pressure distribution with a large aft end loading. This kind of loading 
undercuts the lower surface producing a Whitcomb® type supercritical airfoil. 
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A typical supercritical airfoil design is presented in the Results section. 
The effects of wind tunnel walls in the computation of transonic airfoil 
design and analysis have also been studied. The downstream and upstream 
infinity conditions £or the solldy ideal slotted » and free Jet tunnel walls 
have been derived from a systematic asymptotic solution of the small-distur- 
bance Integrodlfferential equation. 
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SYMBOLS 

scaled mass flux vector 
chord 

lift coefficient 

section normal-force coefficient 
pressure coefficient at half node points 

upper and lower airfoil slopes 
scaled half tunnel wall height 
transonic similarity parameter 
Mach ntimber 

free ctream Mach number 

effective source strength due to airfoil trailing edge openness 

coordinate system 

upper and lower airfoil ordinates 

velocity potential 

far field velocity potential 

angle of attack 

maximum airfoil thickness 


Y 


specific heat ratio 


EQUATION AND BOUNDARY CONDITIONS 

The transonic small disturbance equations are formally derived by an 
asymptotic expansion procedure^ applied to the Euler equations. For an 
airfoil whose upper and lower surfaces are defined by y^ £ ■ 6F^ £(x)-ox the 
perturbation potential satisfies the equation * * 


IK - <Y+1)*.)*„ + ■ 0 


( 1 ) 


The large lateral propagation of transonic disturbances is taken into account 
by the use of the scaled coordinate f ■ y6l/3M2/3, The limit process is 
^ '^hile x,y, and K ■ (1 -m£)/(M^ 62/3) remain fixed. The quantity 

a is the angle of attack. 

Consistent with the small disturbance formulation, the airfoil boundary 
conditions are applied on a slit of 5^ * 0. In the case of pure analysis the 
airfoil boundary condition (flow tangency) is of Neumann type. 


♦y(x,0±) - F^^^(x) - -1 < X 4 1 (2) 


The airfoil leading and trailing edges are at x ■ -1 and x • 1 respectively. 

In the design problem, the airfoil shape corresponding to a given pressure 
distribution Is sought. However, the small disturbance theory cannot resolve 
the nose region accurately. Therefore, the nose shape of an existing airfoil 
is specified for up to 5-lOX of the chord length and a desired pressure distri- 
bution over the rest of the chord Is prescribed. The boundary conditions then 
become a mixed Neumann-Dirlchlet type. On the portion of the airfoil where 
the nose shape Is specified, the boundary condition applied Is given by 
Fq. (2). Over the rest of the airfoil a scaled pressure coefficient 
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- 2 ( 




(3) 


Is prescribed. With known from Eq. (3), the perturbation potential 

^(x,0±) Is calculated by Integration. This value of 4(x,0±) Is then Imposed 
at the airfoil slit as a Dirlchlet type boundary condition. Figure 2 
schematically Illustrates the mixed Neumann-Dirlchlet type boundary condition. 


FAR FIELD 

To avoid mapping procedures which bring Infinity to a finite distance from 
the airfoil, but compromise the difference method, an approximate asymptotic 
solution for ^ valid at large distances from the airfoil Is used as far field 
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boundary condition. The type of far field depends on whether the airfoil is 
kept in free air or a solid, slotted or porous wall wind tunnel. The far 
field expressions to be used in this paper will now be discussed. 


Free Air 


Figure 2 shows the far field arrangement. Along the outer boundary ABODE, 
the perturbation potential 4 Is confuted from 


TF 


n , Qlog? ^ 
2ir ^ 2ir ^ 


(4) 


where V is the circulation around the airfoil, 0 ■ tan“^(»^/x), Q is an 
effect ive sou rce strength due to any openness at the trailing edge, and 
t ■ /x^ + ky^. Only dominant terms are kept in Eq. (4). 


An expression for the source strength Q is obtained by considering the 
integral form of Eq. (1). The divergence theorem is used to obtain this 
Integral relation in the cut region R shown in Figure 3 




(5) 


where B = ^ + tyT is a scaled mass flux vector. Since B is 

is conserved across shock waves, no special boundary terms appear if shocks are 
present in the flow. Expanding Eq. (5) results in 



(rH) 

2 



dy - 0ydx 


( 6 ) 


Substituting Eq. (4) into Eq. (6) and simplifying gives 


- /‘ 

/ic-'-i 


[<|i^]dx, where ( ] denotes Jump 


(7) 
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If the airfoil trailing edge Is closed, then 


/•I 

J [<^^] dx Is zero and no source 


term Is present. Vfhlle designing an airfoil to support a given pressure 
distribution the magnitude of the trailing edge openness Is not knovm a priori. 
In the calculation, Q is therefore evaluated by nonlinear iteration procedure 
analogous to that employed in obtaining the circulation term F for an analysis 
problem. If there Is trailing edge openness, and If the source term Is not 
Included in the far field, then the solution obtained may be questionable. 


Wind Tunnel 


To establish the appropriate boundary conditions on a finite computa- 
tional domain for a tunnel simulation, the far field corresponding to 
X -*• ±» has been derived using a Green's function method. Only the final 
results to the dominant order for the solid, ideal slotted and free jet tunnel 
wall cases are reported here. The airfoil is positioned midway between the 
walls. Thus, the airfoil slit Is at y <■ 0 and the tunnel boundaries are at 
^ ■ ±H. The appropriate boundary conditions on the tunnel walls are 


(x,±H) - 0 

for solid wall 

(x,±H) ± (Kx,±H) » 0 

for slotted 
and free jet 

parameter: 



( 8 ) 


F . ^ ta COS.C (g) 


(9) 


where s Is the distance between slot centers and a is the slot width. For 
free jet case F ■ 0. 

To dominant order in the Karman-Guderley (x,j^) plane, the far fields 
corresponding to these tunnel cases at x ->■ are: 

Solid wall: 

4M - + y\(5)(i? 


2K I 16h^ J 
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where t(l) = ^’^(1) - dummy variables for x and y, h - H, 

H(x) = 1 for X > 0 and H(x) * 0 for x < 0, and u « 4 . 

X 

Slotted and free-jet wall; 


0 


as X -»■ -09 


^ 8gn y - ^ + * * • as X +oa 


Eqs. (11) assume no downwash at upstream infinity. Other expressions can be 
derived for different upstream assumptions. In this connection, the drag. 

lift, and pressure distribution can be shown to be unaffected by addition of 
upstream down flow. ^ 


CONSISTENCY 

When the pressure distribution from the analysis calculation is used as 
an input for design, a consistent discretization procedure would recover the 
airfoil shape exactly, even across the shock. To achieve this agreement, 
du^y points below the airfoil surface are used. Figure 4 shows the grid 
points near the airfoil surface. The points (i, j-1) are the dummy points. 

Analysis 

From the known airfoil shape (♦y)ij,at the beginning of each relaxa- 
tlon cycle the dummy point values <|)i j_i are obtained from the central 
difference formula 




This expression is then used in the finite differenced form of 4, in 
Eq. (1) at the airfoil points (l,j) ?? 


[(k - 
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The nonlinear term In Eq. (13) is central differenced at elliptic points and 
one-sided backward differenced at hyperbolic points. To Improve stability 
near the sonic region Jameson's" pseudo-time operator Is used In the relaxa- 
tion procedure. Once the analysis procedure converges » the pressure 
distribution on the airfoil Is computed at half node points 


C 


P 



‘*'x|l-l/2,J 


( 14 ) 


where 




• Vl.i^ 


( 15 ) 


Design 

As mentioned earlier In this paper, a portion of the airfoil shape 
(<|)y) near the nose is specified (from x - -1 to x “ XD in Figure 2) . On the 

portion of the airfoil under design, the pressure coefficient Cp Is 

i— 1 / 2 , j 

specified at half node points. From Cpj^_i /2 j* perturbation potential 
on the airfoil surface la obtained from Eqs. (14) and (15) 








( 16 ) 


At the airfoil grid points where the shape la specified, Eq. (13) is used 

to evaluate the potential (usual SLOR scheme). Since this potential 

keeps changing during the relaxation cycle, the potential * over the 

accordingly at the beginning of each relaxation 

sloM analysis-design solution converges, the 

slope of the resulting airfoil is computed from Eq. (13) as 




(17) 
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From the slope the airfoil ordinates are 
formula 


* */* 

*D 


calculated from 



the quadrature 


(18) 


Xjj<X<l 


where yu^ji(XD) are the upper and lower y values of the airfoil at X ■ XD, 
where the boundary condition changes from analysis to design. Since Eqs. 
(13), (14), and (15) are used In the same manner In the analysis and design . 
the discretization procedure Is consistent . 


RESULTS 

To check the consistency logic developed in the previous section, a test 
case was run. An analysis calculation for the NACA 0012 airfoil at Moo ■ 0.75, 
Ot ■ 2* was first performed. The pressure distribution from the analysis 
calculation was then used as an Input for the design problem to check if the 
NACA 0012 airfoil shape would be recovered even across the shock. The results 
are shown in Table 1. The upper (yy) and lower (yj^) surface ordinates from 
the design calculation agreed with the original NACA 0012 airfoil up to four 
significant figures, thus establishing consistency of the method. 

Figure 5 shows a design calculation performed on the previous NACA 0012 
analysis solution to get rid of the upper surface shock. The dotted line 
shows the analysis solution, and the solid line shows the prescribed shockless 
pressure distribution. The resulting airfoil Is shown by the solid line. It 
is seen that a slight flattening of the upper surface determined In the design 
phase eliminates the shock. This reshaping produces an openness at the trail- 
ing edge which was properly accounted for by the effective source term In the 
far field. 

Analysis calculations were performed over the shock-free airfoil designed 
In free air, (shown In Figure 5), using solid wall tunnel boundary conditions. 
The airfoil t'ust produces a shockless pressure distribution In free air may 
produce a shock when tested In the wind tunnel. This is Illustrated by the 
results shown In Figure 6. From the figure, It Is evident that when the 
tunnel wall is sufficiently far away from the airfoil, the pressure distribu- 
tion remains shockless. As the tunnel wall Is brought closer to the airfoil, 
the shock appears and moves downstream. 
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Figure 7 shows a supercritical airfoil design with NACA 0012 nose shape. 
First, an analysis solution was generated over the NACA 0012 airfoil at 
Mo# ■ 0.8 and o ■ 2*. This Is shown by the dotted line in Figure 7. Then a 
supercritical pressure distribution was specified with a large loading on the 
aft end of the airfoil. The resulting airfoil resembles a Whitcomb type super- 
critical airfoil which is characterized by a substantially reduced curvature 
on the mid chord region of the upper surface together with increased camber 
near the trailing edge. 

An off design calculation or this supercritical airfoil at " 0.78 
is shown in Figure 8. For qualitative comparison, off design and design 
calculations on the NASA 11% thick supercritical airfoil are also shown. At 
off design Mach numbers the shock reappears. 

The effect of wind tunnel walls on the performance of the free air 
shock free supercritical airfoil is shown in Figure 9. Only the upper surface 
pressure is shown for solid, ideal slotted (F ■ 0.279) and free Jet (F - 0) 
tunnel cases. The solid line in Figure 9 refers to the shock free, free air 
pressure distribution. For tunnel wall height to a chord ratio of 6, the shock 
reappears in all the tunnel wall cases. As expected, the solid wall produces 
a stronger shock while the free jet case produces a weaker one. 


CONCLUDING REMARKS 

An efficient and inexpensive design-analysis code has been developed for 
two-dimensional airfoils. The consistent differencing procedure employed at 
the airfoil boundary allows use of the same code in either analysis or design- 
analysis mode without requiring any modification. In the algorithm, the 
effect of an open trailing edge is properly accounted for, a factor ignored 
by other workers. The fact that no mapping is involved in the two-dimensional 
work makes the extension of the algorithm to three dimensional wing design 
feasible and attractive. 
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Figure 3.- Cut region for far field. 
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Figure 6.- Effect of solid tunnel walls on the performance of 
free-fllght shock-free airfoil. = 0.75; a - 2®. 







Figure 7.- Supercritical-airfoil design with NACA 0012 nose shape. 

M = 0.8; a = 2“. 
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Figure 8.- Off-design pressure distribution. 
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